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The latest experimental data on nuclei at Sn permit 
us for the first time to determine the spin-orbit splittings of 
neutrons and protons in identical orbits in this neutron-rich 
doubly-magic region and compare the case to that of 208 Pb. 
Using the new results, which are now consistent for the two 
neutron-rich doubly magic regions, a theoretical analysis de- 
fines the isotopic dependence of the mean field spin-orbit po- 
tential and leads to a simple explicit expression for the dif- 
ference between the spin-orbit splittings of neutrons and pro- 
tons. The isotopic dependence is explained in the framework 
of different theoretical approaches. 

PACS number(s): 21.60. Cs, 21.10.Pc, 21.60. Jz, 24.10. Jv 



I. INTRODUCTION 

Spin-orbit splitting of the mean field orbitals is one of 
the main factors, which determine nuclear structure in 
nuclei both near and far from the closed shells. While 
the global characteristics of spin-orbit splitting are well 
known, one cannot say the same about the isotopical de- 
pendence of splitting. However, the new experimental 
results obtained recently ]l]-[3) on nuclei close to 132 Sn 
allowed to define a nearly complete set of neutron and 
proton single-particle orbitals and some important stati- 
cal and dynamical properties of the mentioned nuclides. 
In particular, from the measurements on 133 Sb new in- 
formation was obtained Q on the energies and the de- 
cay properties of proton single particle states above the 
Z = 50, N = 82 shells. One purpose of our study is to 
evaluate the new results on the single particle levels at 
132 Sn and intercompare them to the known data on such 
states in the other doubly closed shell (DCS) regions. 
However, our main aim is to examine the magnitude of 
the spin-orbit splittings of neutrons and protons in iden- 
tical orbits, and to determine their isospin dependence. 
A preliminary account of this work was given in p| . 

The determination of the isospin dependence has a 
broader significance, since the magnitude of the spin- 
orbit splittings could be one of the factors contributing 
to significant structural changes in nuclides having an 
extreme neutron excess. Consequently, the trends ex- 
tracted from the empirical data presently available are 
crucial guidelines for the theoretical analysis, and it is 
important to show that model calculations indeed do re- 
produce such trends. 



The following presentation begins with an analysis of 
the existing experimental data in Section II, followed by 
theoretical considerations and evaluations in Section III. 
A summary and conclusions are included in the last sec- 
tion. 



II. EXPERIMENTAL VALUES 

We now examine the experimental data on the spin- 
orbit splitting for neutrons and protons for two groups 
of the doubly closed shell nuclei: the N — Z and the 
neutron-rich regions. In order to facilitate the compari- 
son, the available systematics of single particle energies 
at the N = Z nuclei of 16 O, 40 Ca, 100 Sn, and the neutron- 
rich 132 Sn and 208 Pb are presented in Tables 1 to 5, re- 
spectively. Here, the energies of the particle and hole 
states closest to the Fermi level were determined from 
the differences of binding energies of the core and the 
corresponding adjacent odd nuclei: e (particle) =B(core)- 
B(core+nucleon) and e (hole) =B(core-nucleon)-B (core), 
using the experimental binding energies from [[| . The en- 
ergies of orbitals more distant from the Fermi-level were 
subsequently defined by the addition (subtraction) of the 
experimentally determined excitation energies 
of the corresponding orbitals in the adjacent odd nuclei. 
Important for this process is that we have accounted for 
the fragmentation of states in the cases when the perti- 
nent data were available. The cases where the effect is 
essential are indicated in Tables 1—5 by an asterisk (*) 
next to the experimental value. 



A. N=Z doubly-closed shell regions 

The spin-orbit splittings in the N = Z nuclei can be 
determined using the data in Tables 1-3. Here the spin- 
orbit splittings of the lp\/2— 1?>3/2 an d 1^3/2 — 1^5/2 or- 



bits in 16 O [ |12[]l3|| are practically equal for protons and 
neutrons, with the values of 6.32 and 6.17 MeV (the dif- 
ference being -2.4%) and 5.00 and 5.08 MeV (the dif- 
ference being +1.6%). Similarly, for 40 Ca, Ref. @, we 
have the values for protons and neutrons equal to 2.01 
and 2.00 MeV for the 2p 1/2 -2p 3/2 orbit, 6.00 and 6.00 
MeV for the ld z / 2 -ld b / 2 orbit, and 4.95 and 4.88 MeV 
for the 1/5/2- 1/7/2 case. In the absence of experimental 
data on the single particle states at 100 Sn, we adopt here 
the extrapolated single particle energies from Grawe et 
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al. |l5|, as is shown in Table 3. Based on these data, one 
may conclude that, within the errors, the spin-orbit split- 
ting of the l<77/2— 1<?9/2 orbit is also equal for protons and 
neutrons, namely, at 6.82(28) and 7.00(28) MeV. These 
six cases show that the splitting for the N — Z DCS re- 
gions is practically equal, with small oscillations either 
way, but below 2.6%. This equality simply reflects the 
concept of isobaric invariance in nuclei. 

B. Neutron-rich doubly-closed shell regions 

For the neutron-rich nuclei we find the situation con- 
siderably different, see Tables 4-5 for the data, and Ta- 
ble 6 regarding the splittings. In the 132 Sn region, the 
energy of the 3/2 + proton state in 133 Sb was recently 
determined Q at 2.44 MeV. Using this value and the 
previously determined single particle excitations in nu- 
clei close to 132 Sn (see @J§-§]) the spin-orbit splittings 
of the 2d levels both in proton and neutron systems at 
132 Sn can be now defined. The 2d 3 / 2 — 2d 5 / 2 splitting was 
found to be 1.48 MeV for protons and 1.65 MeV for neu- 
trons. This means that the neutron spin-orbit splitting 
is somewhat larger (by more than 11%) than for protons. 

In the case of 208 Pb, it was noted ||] that the situation 
seemed to be quite opposite. Namely, a simple analysis 
of the single particle levels in 209 Bi and 207 Pb suggested 
Q that the spin-orbit splitting of the 2/ 5 / 2 — 2/ 7 / 2 or- 
bit is equal to 1.93 MeV for protons and 1.77 MeV for 
neutrons. However, a significant correction is needed. It 
follows from the experimental evidence that the neutron 
2/7/2 state in 207 Pb is strongly fragmented, while the 
conclusions in Ref. Q were derived by considering only 
the lowest, albeit the strongest component of this state. 
In order to identify, in the spirit of Refs. |ll|[l7]], the true 
single particle energy of the neutron 2/ 7 / 2 state, we use 
the weighted average of the fragmented 7/2~ energy lev- 
els, with the weight provided by the spectroscopic factors 
determined in the (d,t) reaction on 208 Pb |lq| . In this 
way we obtain a more accurate unfragmented excitation 
energy of this state equal to 2.70 MeV (instead of 2.34 
MeV). Using this excitation energy, included in Table 5, 
we find the neutron spin-orbit splitting of the 2/ orbit as 
2.13 MeV, which, similarly to the case of the 2d orbit in 
132 Sn, is larger by about 10% than the splitting of 1.93 
MeV for protons. 

An additional piece of evidence along the same line is 
given by the analysis of the 3p!/ 2 — 3p 3 / 2 spin-orbit split- 
ting at 208 Pb. One obtains 0.85 MeV for protons (after 
correcting for the fragmentation of the proton 3j?i/ 2 level) 
and 0.90 MeV for neutrons. Thus again the value for 
neutrons is larger by about 6% than for protons. Con- 
sequently, based on the three cases described above, it 
is evident that the neutron spin-orbit splitting in the 
neutron-rich DCS nuclei of 208 Pb and 132 Sn, is system- 
atically larger by ~J0% than the corresponding proton 
splitting. 



C. Fragmentation of strength at Sn 

The fragmentation of single particle states at 208 Pb 
is mainly influenced by the presence of a very low-lying 
and highly-collective 3± phonon state at 2.62 MeV. The 
effects caused by the 2\ state are less as the collec- 
tivization of quadrupole phonon, and the correspond- 
ing nucleon-phonon vertexes are small in heavy nuclei 
near doubly closed shells. For example, in the neutron 
"hole" 207 Pb nuclei, the main part of energy shift of the 
7/2~ level is caused by mixing with the higher lying 
(37/ ® ^l*i3/2 _1 )7/2- state. Numerical evaluation per- 
formed by using the quasiparticle-phonon model with 
the coupling constant extracted from the B(E3; 3^ — > 
ground state) value shows that the 7/2~ level corre- 
sponding to the "pure" v2f- ! / 2 ~ 1 state moves down by 
the amount of ~ 0.4 MeV, thus approaching the ex- 
perimental value of 2.34 MeV. This large shift is due 
to rather strong coupling constant and non spin-flip na- 
ture of the matrix element. The magnitude of the pre- 
dicted shift is very close to the experimental value of 
2.70 - 2.34 = 0.36 MeV mentioned above. At the same 
time, the (37/(g)^lii3/2 _1 ) configuration has no 5/2~ com- 
ponent and thus one does not observe experimentally the 
fragmentation of the lower lying 5/2~ level at 0.57 MeV. 
We note here that the 7/2~ and the 5/2~ states in the 
proton "particle" 209 Bi nucleus have the opposite order- 
ing, 7/2~ being the lower one. Due to mixing with the 
(3^(8)^1^3/2) configuration the "pure" 7r2/ 7 / 2 level is also 
pushed down, but only by about 0.2 MeV due to larger 
energy difference. Thus, after taking account of config- 
uration mixing, not only the neutron A»" (2/) splitting 
between the pure states increased as compared to 1.77 
MeV, but also the proton A^(2/) splitting decreased to 
a smaller value. 

Turning to the region of 132 Sn, we note that the cor- 
responding experimental data on fragmentation of single 
particle states are not known at present. However, as was 
pointed out by Blomqvist |HJ, the 132 Sn and 208 Pb nuclei 
are in some respect twins, having similar shell structures 
with the correspondence of I — ► l + l, j — > j + 1 for most of 
the orbitals in these regions. Therefore, all the arguments 
presented above for splitting of the 2/ levels at 208 Pb are 
completely valid also for the 2d states at 132 Sn, with re- 
placement of l«i3/2 by l/in/2- So far there is no direct ex- 
perimental data on the B(E3; 3|f — > ground state) value 
in 132 Sn. However, the core has much higher rigidity 
here in comparison with 208 Pb and the energy of the 3j~ 
state is substantially higher at 4.35 MeV. Thus from ac- 
counting for configuration mixing one expects some fur- 
ther increase of the (2d) splitting and a decrease of 

A^ (2d), but these changes should be smaller than for 
the 2/ levels at 208 Pb. Estimates based on an indirect 
evaluation of the B(E3) value from the magnitude of the 
octupole effective charge in 134 Te [^0| confirm the pat- 
tern of changes of the A < f s ' n \2d) values presented above. 
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However, in the absence of experimental data on direct 
reactions we present in Tables 4 and 6 the values of en- 
ergies at 132 Sn that do not include averaging over spec- 
troscopic factors. 



III. THEORETICAL APPROACH 
A. General considerations 

Turning to the theoretical interpretation of the ex- 
perimental values of the spin-orbit splitting discussed 
above, we shall first recall that from the point of view of 
many-body theory the average spin-orbit potential has 
its origin in the pair spin-orbit interaction between nu- 
cleons (with tensor forces providing a minor contribution 
as well). On the level of qualitative arguments it was 
noted by Bohr and Mottelson [|l| that due to the symme- 
try properties one should expect the neutron spin-orbit 
splitting to be somewhat larger than that for protons 
in heavier nuclei, simply due to a higher number of like 
particles in the neutron case. However, at that time the 
absence of experimental data did not permit a meaning- 
ful comparison with measurements. With the presently 
available data we can fill this gap, providing also some 
quantitative considerations. 

The two-body spin-orbit interaction differs from zero 
only in the states with a total spin S = 1. The neutron- 
neutron and proton-proton systems have the total isospin 
T = 1 and thus due to the Pauli principle have odd values 
of the relative orbital momentum L (in fact, L = 1). At 
the same time, the neutron-proton system is composed 
from the T = and T = 1 states with equal weights, 
having L = and L = 1, correspondingly. Due to the 
absence of spin-orbit interaction in states with L = 0, 
the pair spin-orbit np interaction is half as strong as that 
in pp or nn-systems. 

If Ut s (n) and Ui s (p) represent the magnitudes of the 
mean spin-orbit field for neutrons and protons and "d(T — 
l,S l =l,L = l)isa quantity representing the parameter 
of the pair spin-orbit interaction in a state with T = 
1, S = 1, L = 1 then the above discourse leads to 



U ta {n) ~ 0(1, 1, 1) • NV + -Z = • [A - 
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and 



U u (p) ~ 0(1, 1A)-(j+z\=^-(a-j\. (1) 

As the spin-orbit splitting A^™' p ^ ~ Ui s [n,p), the rel- 
ative difference "e" of the neutron and proton spin-orbit 
splittings is given by the expression: 
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On the other hand, we can express the strength of the 
spin-orbit mean field in the form: 



1 N — 

U es (r 3 ) = V es [\ + -p u — 



(3) 



Here T3 = — 1 for neutrons, r 3 = +1 for protons and Pi s is 
the parameter that defines the isospin dependence of the 
mean spin-orbit field. Then we easily obtain, this time in 
terms of eq. (3), an expression for the relative difference 
between the spin-orbit splittings of neutrons and protons 
in identical orbits, e: 
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It follows from a comparison of eqs. (2) and (4) that 
Pu = -2/3. 

Strictly speaking, this derivation was performed for the 
two-body spin-orbit interaction. However, as mentioned 
above, tensor forces provide also some contribution to 
the spin-orbit splitting. This non-central interaction is 
proportional to S12 with 



S12 = 3((Tin)((T2n) - eri<7 2 



[[(71 ® <T2? <8> Y 2 }° . 



(5) 



One can easily see from (5) that the diagonal matrix el- 
ements of this interaction are different from zero only for 
states with S — 1 and L > 1 , of which the S = T = L = 1 
one is of the main importance. It is just the state which 
was already considered in this subsection in the case of 
spin-orbit interaction. Consequently, the diagonal part of 
tensor forces also provides contribution of the type given 
by eq. (3) with j3g s = —2/3, and thus it leads only to a 
renormalization of the Vi s value. However, as the spa- 
tial part of tensor operator is proportional to ^(n) and 
due to the spin structure of S12, this renormalization 
equals zero in cases of spin saturated spherical nuclei. 
Thus in 16 O and 40 Ca tensor forces give a contribution 
to the isoscalar part of the spin-orbit splitting, that is 
mediated by their non-diagonal part and caused by ad- 
mixtures, that are out of the Hartree-Fock type ground 
state. As was shown in Ref. [ p2| , tensor forces may re- 
ally lead to a substantial contribution to the isoscalar 
part of spin-orbit splitting. At the same time, in nuclei 
that are not spin saturated, such as 48 Ca, tensor forces 
can contribute to the spin-orbit splitting even in the "di- 
agonal" scheme (i.e.: a scheme without admixtures), if 
the antisymmetrization is properly included. Our nu- 
merical calculations for seniority one states of 47 Ca and 
47 K both having one neutron or proton hole and per- 
formed in the framework of the multiparticlc shell model 
with tensor forces taken from our previous works p3|]- 
[E7fl, have demonstrated that the inclusion of a tensor 
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component of the interaction leads to energy shifts that 
correspond to some variation of the spin-orbit splittings 
A ta , such that in 48 Ca A^(ld) - A^(ld) = 0.34 MeV 

and A^ } (lp) - A^(lp) = 0.24 MeV. These shifts arise 
from neutrons filling the vlf-j/2 subshell and are mainly 
due to charge exchange two-body matrix elements of the 
np-interaction mediated by the isovector part of the ten- 
sor force (~ T1T2). Thus the inclusion of tensor forces 
does not change the pattern of spin-orbit splitting, which 
also leads to negative values of (3i s ranging from about 
—0.4 to —0.7. These results qualitatively agree with 
those presented in Ref. pS[ ], where in the framework of 
the Brueckner-Hartree-Fock method with Reid potential 
(containing both the spin-orbit and tensor components), 
a substantially larger neutron than proton splitting was 
obtained for the lp and Id orbitals in 48 Ca with /3i s in 
the range from about —0.5 to —1.8. We note that the 
data on spin-orbit splittings of the 2d states in 132 Sn as 
well as on the splittings of the 2/ and 3p levels in 208 Pb 
lead to effective values of fy s equal to —0.55, —0.60 and 
—0.27, respectively, which are numbers in very satisfac- 
tory general agreement with the prediction of eq. (4). 

It is thus of substantial interest to evaluate to what 
extent the isotopic dependence of the spin-orbit splittings 
are reproduced by standard model calculations. Three 
different approaches were made as described below. 



B. Evaluation I: Walecka model 

The first evaluation is made in the Hartree approxima- 
tion starting from the Dirac phenomenology with meson- 
nucleon interactions according to the Walecka model [0 . 
One obtains (see for example 1 30 1 — |57j and references 
therein) a Skyrme-type single particle equation for a nu- 
cleon having the effective mass m* N . This approach well 
explains the magnitude of spin-orbit splitting in nuclei. 
Here, and mainly for heavier nuclei, we concentrate only 
on the difference between the proton and neutron split- 
tings of spin-orbit partners in the same nuclei and result- 
ing from a spin-orbit potential having the form (see for 
example p§- pi): 



77 — ^ N 1 f(—H\ 2 ^ m/O Q° \ 



where V = V° - r 3 • V 1 and S = S° - t 3 • S 1 are the 
vector and scalar fields related to corresponding mesons, 
m* N = mN+\(S—V), while k is the ratio of tensor to vec- 
tor coupling constants of p-meson. Various approaches 
have been used to determine the coupling constants. In 
p6| the meson-nucleon coupling constants, defining the 



V and S fields, were taken from the Bonn NN boson ex- 
change potential fl38|| , where <r and tro are scalar mesons 
imitating the 2ir exchange in the NN- systems with T=l 
and T=0, correspondingly. In other works (see for exam- 
ple J33|- p5| ) the constants were defined from the de- 
scription of global nuclear properties, with inclusion of 
the <7 3 and o terms in the Lagrangian density (one <r- 
meson with the same characteristics for T=l and T=0 
channels was used, which leads to zero contribution of 
this meson to S 1 in formula (6); note also that the ten- 
sor term was not included in the p-meson vertex in Refs. 

HI" II); 

Taking into account that the radial dependence of the 
(jbjv/jtijv) is much weaker than that of V and 5, which 
are considered to be proportional to the density in the 
form of Fermi function, one can approximately present 
formula (6) as follows: 



\_df_ 

r dr 



N ■ 



/ = [1 + exp(- 



A 



R 



■ 7 3 



£ -s: 



)]■ 



(7) 



Calculating the V and S magnitudes in the center of 
nuclei at the values of vector and scalar densities p v — 
0.17, p s = 0.16, p~ = 0.17 (V - Z)/A, pj = 0.16 {N - 
Z) I A (all in fm~ 3 ) , using the coupling parameters from 
p6| , [ p8[ and taking into account the isotopic dependence 
of m^jm* Nl we obtain Vg s « 34 MeV- fm 2 and fii s » 

- 0.40. If we use the NL2 set of parameters from j$4 35 
then we have V ls w 31 MeV • fm 2 , ts w - 0.43. At the 
same time the set NL1 from [ p3p5| ] , giving small values of 
effective masses, leads to Ve 8 ~ 50 MeV • fm 2 and (its ~ 

- 1.3. As the V j S* 1 magnitudes are proportional to p~ 
and pj, both the formulae (6) and (7) give equal spin- 
orbit splitting for protons and neutrons in the N = Z 
nuclei. It should be noted, that the value of (if s is always 
negative and is determined mainly, or entirely, by the 
p-meson conribution. 

The magnitudes of the empirical effective values of (3e s 
at 132 Sn and 208 Pb, listed in subsection IIIA, are quite 
well reproduced by the model calculations in this sub- 
section, in particular by those using the first two sets of 
parameters. 

It is worth mentioning that a study of the neutron 
spin-orbit splitting in light nuclei as a function of A at 
given Z was recently performed in the framework of the 
Walecka model by Lalazissis et al. [[59). However, the 
intercomparison between the splittings of both proton 
and neutron "similar" spin-orbit doublets in the same 
nuclei was not performed there. 



C. Evaluation II: Woods-Saxon model 

In the second approach, using a Woods-Saxon (W-S) 
model, we let the single particle levels be generated by 
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the potential 



U(r,a,T 3 ) = U (T 3 )f(r) + 



Uu(n) df i t 

r dr 



(i + n) 



u, 



Coul i 



(8) 



where U (t 3 ) = V (l + ±/3^ • r 3 ); U ls and f(r,a,R) 
are defined by eqs. (3) and (7), R = r^A 1 ^, while 
Ucoui{r, R c , Z) represents the potential of a uniformly 
charged sphere with the charge Z and radius R c — 
rcA 1 / 3 . 

In previous works |2j| - [^7|, calculations were made 
using the V = -51.5 MeV, r = 1.27 fm, V ls =33.2 
MeV • fm 2 , a(p)= 0.67 fm, a(n) = 0.55 fm and f3 ls =f3 
=1.39, which on the average described the spectra of sin- 
gle particle states in nuclei from 16 to 208 Pb. This set 
of parameters is denoted here as the "Standard" one. 
With the appearance of new experimental data on the 
single-particle levels, we performed a new determination 
of parameter values through the Nelder-Mead method 
by minimizing the root-mean square deviation 
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The computation demonstrated a very small sensitivity 
of results to the value of r c , which was adopted to be 
the same as before: r c = 1.25 fm. The minimization of 
S performed for all nuclei presented in Tables 1-5 with 
r c = 1.25 fm and different values of ro, showed that the 
minimum in all cases corresponds to ro ~ 1.27 fm, that 
also coincides with the value adopted by us before. The 
values r c = 1.25 fm and ro = 1.27 fm were thus fixed in 
further calculations. 

As was noted above, the optimal relation of proton to 
neutron spin-orbit splitting corresponds to fli s ~ —0.6. 
The fourth column, "Set 1", of Tables 4 and 5 presents 
the values of theoretical energy levels obtained in the 
optimization with fixed values of Pi s = —0.6, a p = 0.67 
fm and a n = 0.55 fm. 

The fifth column, "Set 2", of Tables 4 and 5 presents 
the results of optimization with only two fixed parame- 
ters: a p = 0.67 fm and a n = 0.55 fm. 

The values of " Set 3" corresponds to an optimization 
at fixed (3i s = —0.6, while "Set 4" are the results with no 
parameters fixed. 

We see that the optimized values of Vb, Vt s and (3 (see 
formula (8)) are very close to the "Standard" ones, with 
small variations from nucleus to nucleus. The magni- 
tudes of the diffusinesses " a" vary more strongly, differ- 
ing by about 10% to 15% from their "standard" values. 
A comparison of the "Stnd" with "Set 1" and of "Set 
3" with "Set 4" results shows that the contribution of 
(3e s to the root-mean square deviation <5 is small. It is 
thus more reasonable to define (3^ s not from a minimiza- 
tion of 5, but rather by using the experimental and the- 
oretical arguments mentioned above. This conclusion is 



confirmed by the results of Koura and Yamada []4l|] , who 
made a number of different fits of W-S parameters to 
the same set of experimental data, obtaining diverse (in 
magnitude and sign) values of the parameter that defines 
the contribution to the spin-orbit term, which is linear 
in (N — Z)/A. A global adjustment of W-S parameters 
simply appears to be only weakly sensitive to details of 
the spin-orbit splitting. 

As mentioned previously, the energies of levels in nuclei 
with N = Z (see Tables 1-3) are independent of (3 and 
(3i s . Here the optimization was performed twice, first 
with fixed values of a n = 0.55 fm and a p — 0.67 fm 
with a subsequent definition of V and Vi s ("Set 1") and 
secondly without fixing any parameters (" Set 3" ) . 

The results of the calculations presented in Tables 1 
to 5 include some levels having positive energies, i.e. un- 
bound but sub-barrier states. In such cases we present 
here the real part of the single particle energies only for 
those states having very small decay widths. 

To summarize Evaluation II, we have determined the 
parameters of the W-S potential using a global mean 
square-root optimization, except for the isospin depen- 
dent spin-orbit term, where the parameter value was 
found to be insensitive to the adjustment. Hence the 
value of (3e s ~— 0.6 was deduced from physical considera- 
tions based on experimental spin-orbit splittings. 



D. Evaluation III: Hartree-Fock with a Skyrme 
interaction 

For the third model approach, which complements the 
first two evaluations using the empirically adjusted W-S 
potential (8) and the microscopical procedure, we have 
selected the Hartree-Fock calculations with the SIII inter- 
action. The results of these self-consistent calculations, 
listed in the last two columns of Tables 1 to 5, were ob- 
tained by considering the contribution of a single-particle 
part of the center-of-mass energy and taking into account 
the Coulomb exchange term in the Slater approximation. 
The SIII-1 results correspond to calculations which take 
into account all terms of the energy functional contribut- 
ing to spin-orbit splitting, while the SIII-2 results have 
been obtained by omitting the spin density terms in the 
spin-orbit potential. In the last case our results are close 
to those from the study by Leander et al. p2| performed 
for 208 Pb, 132 Sn and 100 Sn nuclei. We see that the results 
obtained in the framework of the Hartree-Fock method 
also demonstrate that the calculated neutron spin-orbit 
splittings of the 2d orbit in 132 Sn as well as of the 2/ 
and 3p orbits in 208 Pb are larger than for protons and 
they correspond to effective j3n 8 in the interval of —0.9 to 
—0.6. We note that the difference between the neutron 
and proton spin-orbit splittings is reproduced here by 
using a simple parameterization of Skyrme forces. Our 
calculated results differ from those of Noble |n| who pro- 
posed that the isotopic dependence of the spin-orbit po- 
tential in the Hartree scheme is cancelled through the 
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contribution of exchange terms, but agree with that of 
|p8| . We mention here that the SIII parameterization 
contains density-dependent terms that imitate in some 
sense the three-body interaction. 

Tables 1 and 2 give results for 16 O and 40 Ca, which 
are spin-saturated nuclei. In these cases the spin den- 
sity terms, included in SIII-1 but not in SIII-2, do not 
contribute significantly to the spin-orbit splitting (the 
contributions in these cases are only due to small differ- 
ences in the radial wave functions of spin-orbit partners) . 
Consequently, as can be expected, the SIII-1 and SIII-2 
calculations give very similar results in both cases. 



IV. DISCUSSION AND CONCLUSIONS 

Using theoretical analysis and systematics of available 
experimental data we have derived formula (3) that de- 
scribes the difference between the neutron and proton 
spin-orbit splittings, i.e. the isotopic dependence of the 
mean field spin-orbit splitting. The splitting becomes 
larger for neutrons than for protons in nuclei having 
N > Z. The general arguments presented initially (based 
on the properties of the two-body spin-orbit and tensor 
interactions) gave a result in fair agreement with the em- 
pirical observations. A further microscopic study within 
the Walecka model supports this initial result, while it 
was found that a global fit of Woods-Saxon model pa- 
rameters appears to be rather insensitive to the isotopic 
dependence of the spin-orbit splitting. A self-consistent 
calculation using the SIII interaction gave results in gen- 
eral agreement with the experiment and prediction by 
eq. (3) with negative values of fy s ~ —0.6. 

In this context, one should point out that within the 
Walecka model, the sign of the isospin term in the spin- 
orbit potential is in agreement with the sign of an analo- 
gous term present in the expression for the central nuclear 
potential. While the spin-orbit term in this model is de- 
fined, very approximately, by the (V — S) combination of 
the entering fields, the central nuclear potential is propor- 
tional to the (V + S) combination. The main, isoscalar, 
part of the (V — S) term is positive and the addition of 
an isovector contribution, arising from V 1 , leads for the 
N > Z nuclei, as was shown above by us, to an additional 
term (positive for neutrons and negative for protons), its 
magnitude growing with (N— Z), together with the ratio 
of neutron to proton splittings. At the same time, the 
isoscalar part of the central (V + S) term is negative. The 
addition of a V 1 term leads here for neutrons in (N > Z) 
nuclei to reduction of the absolute value of (V + S). So, 
with increasing A at a given Z, the depth of the central 
nuclear potential for neutrons decreases and they become 
less bound, while the protons become more bound. All 
this is reflected in the W-S model (see eq. (8) above) by 
the fact that j3i s is negative, while j3 is positive. The two 
models are thus fully consistent in this respect. 

The isotopic dependence of the spin-orbit splitting has 



also been studied with methods somewhat different than 
those used here. In the work of Mairle j43j the average 
spin-orbit potential was obtained as a convolution with 
proton and neutron densities taken in the ratio defined 
by the short-range two-body spin-orbit interaction. How- 
ever the isotopic dependence of the average spin-orbit 
potential was not derived here in an explicit form. This 
point has some importance, since our analysis, based on 
the existing empirical data and different theoretical ap- 
proaches, resulting in a simple expression, immediately 
shows that the difference between the neutron and pro- 
ton splittings becomes saturated at large A, which pre- 
cludes very large differences. The rather modest differ- 
ence with a magnitude of about 10% seen in the 132 Sn 
region is already about 25% of the saturation value, sug- 
gesting that the isospin dependence in itself is unlikely to 
lead to dramatic structural changes. However, in cases 
of extreme neutron excess, when the difference between 
neutron and proton spin-orbit splittings approaches the 
maximum value of about 40% (corresponding to several 
hundreds of keV) a rather significant effect on the order- 
ing of levels can be expected. 
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Table 1. Single particle levels of ib O 



Poschl, P.Ring, 
308 (1967). 



n£j 



Stnd Set 1 Set 3 SIII-1 SIII-2 



'Id- 



v2s l 



3/2 



'Id 



/2 



5/2 
^lP3/2 



(0.94) 
-3.27 
-4.14 
-15.67 
(-21.84) 



0.89 

-3.59 

-6.97 

15.06 

■19.98 



0.18 

-3.89 

-6.85 

16.05 

20.25 



0.20 

-3.31 

-6.41 

16.33 

■20.10 



0.66 

-2.88 

-6.87 

14.58 

20.58 



0.67 

-2.87 

-6.89 

14.56 

20.59 



71-14/2 
7t2si/ 2 

7rlrf 5 / 2 

1"lPl/2 
7TlP3/2 



(4.40) 
-0.11 
-0.60 
-12.13 
(-18.45) 



3.76 
-0.89 
-2.76 
-9.95 
■14.66 



2.92 

-1.14 

-2.67 

■10.87 

14.90 



3.48 

0.22 

-2.97 

12.60 

16.40 



3.55 

0.03 

-3.57 

•11.17 

17.07 



3.56 

0.03 

-3.59 

11.15 

17.08 



Set 1: V = -52.21 MeV, V ls = 28.6 MeV ■ fm 2 , a p = 0.67 fm, 
a n = 0.55 fm are fixed. 

Set 3: V = -51.40 MeV, V ls = 25.7 MeV ■ fm 2 , a p = 0.45 fm, 
a n = 0.50 fm. 



Table 2. Single particle states of 40 Ca. 



n£j 



Stnd Set 1 Set 3 SIII-1 SIII-2 



fl/5/2 

fl/7/2 



vld 



v2s lt 



3/2 



-3.48 

-4.42 

-6.42 

-8.36 

-15.64 

-18.11 

■21.64* 



-2.57 

-3.35 

-5.71 

-10.43 

-16.21 

-16.51 

-21.08 



-3.91 

-4.08 

-6.08 

10.44 

17.40 

17.17 

21.44 



-3.54 

-4.69 

-6.57 

-9.72 

16.43 

17.00 

20.52 



-1.49 
-2.20 
-4.09 
-9.92 
15.53 
15.94 
21.90 



-1.48 

-2.23 

-4.05 

-9.94 

-15.54 

-15.92 

-21.90 



Tl/5/2 

7r2p 3 /2 
Tl/7/2 

7Tlrf 3/2 

Trld 5/2 



3.86 

2.64 

0.63 

-1.09 

-8.33 

-10.85 

14.33* 



4.92 
2.62 
0.89 
-2.19 
-7.11 
-8.18 
-12.05 



3.79 
2.11 
0.60 
-2.18 
-8.25 
-8.78 
12.36 



3.41 
2.07 
0.45 
-2.85 
-9.01 
-9.30 
13.19 



4.90 
3.66 
2.23 
-3.04 
-8.52 
-8.77 
14.74 



4.91 
3.64 
2.26 
-3.06 
-8.53 
-8.75 
14.75 



"Set 1": V = -52.39 MeV, V u =27.9 MeV ■ fm 2 ; a p = 0.67 fm 
and a n = 0.55 fm are fixed. 

"Set 3": V = -52.95 MeV, V ls =28.2 MeV ■ fm 2 , a p = 0.63 fm, 
a n = 0.68 fm. 

Experimental single particle energy marked by an asterisk (*) rep- 
resents a mean value weighted by the spectroscopic factors. 
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Table 3. Single particle states of luu Sn 



Table 5. Single particle states of 20 «Pb 



nij 



Stnd Set 1 Set 3 SIII-1 SIII-2 



nlj 



^exp 



Stnd Set 1 Set 2 Set 3 Set 4 SIII-1 SIII-2 



vlhi 
v2d- 



1/2 
3/2 

"1-97/2 

v2d 5 / 2 



-8.6(5) 
-9.2(5) 
-9.3(5) 
-10.93(20) 
-11.13(20) 
-17.93(20) 
-18.38(20) 



-8.66 
-8.90 
-9.16 
11.64 
11.62 
17.23 
19.14 



-9.01 

-9.24 

-9.53 

12.02 

11.97 

17.61 

19.53 



-8.72 
-8.70 
-9.13 
11.23 
11.59 
17.21 
18.93 



-6.35 

-7.84 

-7.58 

10.33 

10.07 

16.54 

19.08 



-6.87 

-7.66 

-7.52 

-9.63 

10.10 

17.00 

18.93 



7rlg 7/2 
7r2d 5/2 

1159/2 
7l2}33/2 

711/5/2 



3.90(15) 


3.88 


3.54 


2.70 


3.38 


4.04 


3.00(80) 


2.74 


2.45 


2.64 


3.70 


3.69 


-2.92(20) 


-2.01 


-2.36 


-3.66 


-2.74 


-3.16 


-3.53(20) 


-3.48 


-3.84 


-3.94 


-4.80 


-4.65 


-6.38 


-4.95 


-5.31 


-5.55 


-6.22 


-6.18 


-8.71 


-5.54 


-5.92 


-7.60 


-8.43 


-7.89 



Set 1: Vb = -51.97 MeV, V ls =33.5 MeV • fm 2 ; a p = 0.67 fm and 
a n = 0.55 fm are fixed. 

Set 3: Vb = -51.40 MeV, V la =35.6 MeV ■ fm 2 , a p = 0.52 fm, 
a n = 0.56 fm. 



Table 4. Single particle states of 16 Sn 



Stnd Set 1 Set 2 Set 3 Set 4 SIII-1 SIII-2 



^2/5/2 

ulhg/2 

^2/7/2 
i/2d 3/ 2 

^1^11/2 
u3s 1/2 

v2d 5/2 



-0.58 
(-0.92) 
-1.02 
-1.73 
-2.58 
-7.31 
-7.55 
-7.64 
-8.96 
-9.74 



0.36 
-0.13 
-1.61 
-0.78 
-2.18 
-7.74 
-7.11 
-7.68 
-9.66 
-10.56 



0.73 
-0.48 
-0.84 
-0.88 
-2.55 
-7.45 
-7.96 
-7.73 
-9.94 
10.04 



0.46 
-0.09 
-1.38 
-0.77 
-2.21 
-7.62 
-7.23 
-7.64 
-9.66 
10.39 



0.22 
-0.55 
-0.47 
-1.42 
-2.84 
-7.63 
-7.33 
-8.03 
-9.98 
-9.51 



-0.01 
-0.61 
-0.97 
-1.32 
-2.52 
-7.77 
-6.60 
-7.93 
-9.69 
-9.81 



0.67 

0.16 

-0.72 

-0.16 

-1.67 

-8.42 

-7.69 

-8.26 

-10.71 

-11.92 



0.79 

0.20 

-0.02 

-0.14 

-1.71 

-8.26 

-8.23 

-8.21 

10.71 

11.32 



n3s 
irlh 
7r2d : 
7T2d. 



1/2 
11/2 



3/2 
5/2 



T197/2 

T199/2 
7r2p 1/2 



(-6.83) 
-6.84 
-7.19 
-8.67 
-9.63 
-15.71 
-16.07 



-6.84 
-7.32 
-6.86 
-9.36 
-9.84 
-14.91 
■16.01 



-6.87 
-6.66 
-7.20 
-9.20 
10.41 
-14.46 
16.22 



-6.80 
-7.46 
-6.74 
-9.37 
-9.66 
15.00 
15.92 



-6.64 
-6.77 
-7.07 
-9.04 
10.60 
14.57 
16.14 



-6.70 
-7.48 
-6.72 
-9.30 
-9.81 
15.02 
15.91 



-4.97 
-5.64 
-5.93 
-7.88 
10.08 
15.03 
16.68 



-4.90 
-6.01 
-5.77 
-7.88 
-9.56 
15.36 
16.55 



"Stnd": 8 = 0.589 MeV. 

"Set 1": Vb = -51.56 MeV, V ts 

8 = 0.638 MeV. 

"Set 2": V = -51.44 MeV, V ta 
f3 ts = 1.35, 8 = 0.575 MeV. 
"Set 3": Vb = -51.55 MeV, 



33.3 MeV ■ fm 2 , 



34.8 MeV ■ fm 2 



fm 2 



a p = 0.63 fm, a„ 
"Set 4": Vb = 



1.39, 



1.39, 



1.31, 



1.34, 



V ts = 32.4 MeV 
= 0.66 fm, 8 = 0.546 MeV. 
-51.56 MeV, V la = 34.1 MeV ■ fm 2 
P la = 1.33, a p = 0.65 fm, a n = 0.66 fm, 8 = 0.478 MeV. 
Note that some theoretical works [Q] postulate that the neutron 
1*13/2 state at 132 Sn is only 1.9 MeV above the ^2/ 7 / 2 level. Our 
calculations unequivocally demonstrate, that this state lies consid- 
erably higher, with it's energy equal to +0.55, +1.59 and +1.02 
MeV for the "Stnd", SIII-1 and SIII-2 parameter sets, respectively. 



u3d 



'3/2 

^2g 7 /2 



u4s 



1/2 



vZd 



vXi 



5/2 



11/2 
^2g 9/2 
i/3p 1/2 
^2/ 5 /2 



vXi 



13/2 
^2/7/2 



ulh 



9/2 



-1.40 
-1.44 
-1.90 
-2.09* 
-2.37 
-3.16 
-3.94 
-7.37 
-7.94 
-8.27 
-9.00 
-10.07* 
-10.78 



-0.32 
-0.79 
-0.80 
-2.42 
-1.50 
-4.24 
-3.71 
-7.32 
-8.42 
-8.18 
-9.21 
-10.57 
-12.06 



-0.02 
-0.18 
-0.70 
-3.05 
-1.45 
-3.37 
-3.82 
-6.94 
-7.87 
-8.03 
-9.62 
10.57 
11.35 



-0.23 
-0.65 
-0.74 
-2.31 
-1.40 
-4.05 
-3.59 
-7.17 
-8.25 
-8.04 
-9.08 
10.43 
11.87 



-0.96 
-0.89 
-1.63 
-2.23 
-2.35 
-2.71 
-4.24 
-7.59 
-8.17 
-8.59 
-8.84 
10.72 
10.60 



-0.99 
-1.14 
-1.51 
-1.55 
-2.13 
-3.33 
-3.88 
-7.61 
-8.38 
-8.43 
-8.31 
10.46 
11.09 



0.38 

0.01 

-0.08 

-1.41 

-0.39 

-3.37 

-2.91 

-7.21 

-8.59 

-8.18 

-9.73 

11.21 

13.16 



0.42 

0.14 

-0.06 

-1.93 

-0.38 

-2.77 

-2.97 

-7.13 

-8.44 

-8.15 

10.21 

11.24 

12.67 



7r3p 1/2 
T3p 3/ 2 
""2/5/2 



13/2 
t2 /7/2 



wlh. 



9/2 
7r3s 1/2 



-2d- 



irlh 



3/2 



-2d 



11/2 



'5/2 
rl97/2 



0.17* 
-0.68 
-0.97 
-2.19 
-2.90 
-3.80 
-8.01 
-8.36 
-9.36 
-10.04* 
-12.18* 



0.63 

-0.45 

-0.68 

-2.86 

-3.38 

-4.60 

-7.76 

-8.41 

-9.33 

-10.10 

-12.07 



0.43 

-0.46 

-1.03 

-2.37 

-3.24 

-5.11 

-7.86 

-8.66 

-8.99 

10.05 

12.45 



0.72 
-0.35 
-0.60 
-2.71 
-3.26 
-4.53 
-7.67 
-8.32 
-9.18 
-9.98 
-11.99 



0.29 
-0.58 
-1.03 
-1.94 
-3.21 
-4.71 
-7.87 
-8.59 
-8.60 
-9.96 
12.08 



0.47 

-0.69 

-0.61 

-2.78 

-3.53 

-4.01 

-7.87 

-8.30 

-9.21 

10.15 

11.58 



2.79 

1.99 

0.60 

-1.20 

-1.64 

-4.68 

-7.39 

-8.64 

-9.35 

10.29 

13.94 



2.88 

2.03 

0.74 

-1.53 

-1.66 

-4.24 

-7.33 

-8.51 

-9.65 

10.28 

13.59 



The "standard" set of parameters corresponds to Vb 



= -51.50 
0.67 fm, 



MeV, Vis = 33.2 MeV ■ fm 2 , /3 = ft> s = +1.39, a p 
a n = 0.55 fm and 8 = 0.604 MeV. 

"Set 1" corresponds to V = -51.39 MeV, V la = 33.1 MeV ■ fm 2 , 
P = 1.43 with (3 ls = -0.6, a p = 0.67 fm, o„ = 0.55 fm fixed; 
8 = 0.654 MeV. 

"Set 2" corresponds to V = -51.34 MeV, V ia = 33.1 MeV ■ fm 2 , 



1.26 with a„ 



(3 = 1.40, A_ 
8 = 0.593 MeV. 
"Set 3" corresponds to Vb 
13 = 1.36, a p = 0.73 fm, 
(3( s = —0.6 is fixed. 
"Set 4" corresponds to Vq 



0.67 fm, a n 



0.55 fm fixed; 



= -51.99 MeV, V ls = 32.7 MeV ■ fm 2 , 
a„ = 0.72 fm with <5 = 0.369 MeV; 

= -51.93 MeV, V ls = 35.2 MeV ■ fm 2 , 
0.73 fm, a n = 0.72 fm; <5 = 0.366 MeV. 



(3 = 1.38, (3 ls = 1.76, a p 
Experimental single particle energy marked by an asterisk (*) rep- 
resents a mean value weighted by the spectroscopic factors. 

Table 6. Magnitudes in MeV of neutron and 
proton spin-orbit splittings. 

Nucleus nij A exp Stnd Set 1 Set 3 SIII-1 SlTi-2~ 



Sn 



v2d 1.65 
7r2d 1.48 



1.92 
2.50 



2.49 
2.00 



21 )K 



Pb 



v2f 2.13 2.15 2.70 
tt2/ 1.93 2.70 2.21 



vip 0.90 0.86 
n3p 0.85 1.08 



1.09 
0.89 



2.35 
1.97 

2.55 
2.18 

1.00 
0.87 



2.29 
1.95 

2.62 
2.24 

0.97 
0.80 



2.45 
2.11 

2.80 
2.40 

1.02 
0.85 



Notation is as in previous tables. Data are given only in the cases 
where spin-orbit partners of both neutrons and protons in identical 
orbits have been observed experimentally. Note that the splittings 
are practically identical for neutrons and protons in the N = Z 
nuclei, which are not included in this Table. 
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